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Abstract 

We address the problem of Gribov copies in lattice QCD. The gluon propagator is computed, in the 
Landau gauge, using 302 {f3 — 5.8) 12* configurations gauge fixed to different copies. The results of the 
simulation shows that: i) the effect of Gribov copies is small (less than 10%); ii) Gribov copies change 
essentially the lowest momenta components (g < 2.6 GeV); iii) within the statistical accuracy of our 
simulation, the effect of Gribov copies is resolved if statistical errors are multiplied by a factor of two or 
three. Moreover, when modelling the gluon propagator, different sets of Gribov copies produce different 
sets of parameters not, necessarily, compatible within one standard deviation. Finally, our data supports 
a gluon propagator which, for large momenta, behaves like a massive gluon propagator with a mass of 
1.1 GeV. 
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1 Introduction and Motivation 



Quantum Chroniodynamics (QCD) is the theory that describes the interaction between quarks and gluons. 
The definition of the QCD generating functional d la Faddeev-Popov m3 |S| requires a choice of a gauge 
condition, uniquely satisfied in each gauge orbit, i.e. on each set of fields related by a gauge transformation. 
For the Landau, the Coulomb gauge and for small field amplitudes, the gauge condition is uniquely satisfied 
in each gauge orbit. However, if large field amplitudes are involved, the gauge fixing condition has multiple 
solutions in each gauge orbit ^ |S] , the Gribov copies. 

Gribov copies appear when large field amplitudes are involved and rise the question of how to define 
the generating functional for the nonperturbative regime of Quantum Chromodynamics. Moreover, in |^ it 
was proved that it is not possible to find a local continuous and unambigous gauge fixing condition for any 
SU(N) gauge theory defined on the manifold S4. A similar result for the four-torus was obtained in [J]. 

*email: psilva@teor.fis.uc.pt 
temail; orIando@teor.fis.uc.pt 



1 



For the continuum formulation of QCD, Zwanziger argued in |S] that the Landau gauge Faddeev-Popov 
formula 

5{dA) dei[-d-D{A)] eM-SYM{A)] , (1) 

restricted to the region where the Faddeev-Popov operator is positive definite —d ■ D{A) > 0, the Gribov 
region fi, provides an exact non-perturbative quantisation for QCD. This result helps to eliminate some 
theoretical questions about the investigations of QCD using Dyson-Schwinger equations (DSE). Nevertheless, 
in what concerns the non-pcrturbative regime of QCD, being unable to solve exactly the DSE, the results of 
such studies should be compared to lattice results. In this way one can test the validity of the approximations 
and ansatz used to solve the DSE and, simultaneously, the lattice algorithms. 

The formulation of gauge theories on the lattice does not require gauge fixing. As long as one is interested 
only on gauge invariant operators, the lattice calculation is not plagued with the problem of Gribov copies. 
However, the investigation of the Green's functions of the fundamental fields, such as the gluon, ghost and 
fermion propagators, implies the choice of a gauge. On the lattice, typically, a simulation begins by generating 
a number of thermalized gauge configurations. In order to compute, for example, the propagators, each 
configuration is then rotated to satisfy a given gauge fixing condition. Finally, the propagator is computed 
using these rotated configurations. For the Landau gauge, gauge fixing is implemented by computing a 
maximum of a given function defined on the gauge orbits. Now, the problem of the Gribov copies is due to 
the several maxima of the maximizing function. The first observations and studies of lattice Gribov copies 
were done long ago |U1 IIUL [TTl 1121 113L I14| . However, how the choice of Gribov copies change the correlation 
functions is not yet clear. 

On the continuum formulation, Gribov 4 studied SU(2) gauge theory. His purposal to solve the 
problem of the different copies was to restrict the functional integration space to the so called Gribov region 
ri. The gluon propagator computed by functional integrating the gluon fields over f2 does not show the usual 
perturbative l/q^ behaviour but, instead, /{q^ + M^), with M being a mass scale which measures the 
volume of Vl. Note that the two propagators agree for the high energy regime. 

On the lattice, there was a number of studies about Gribov copies and different observables in various 
gauges. In this paper we will be mainly concerned about the gluon propagator computed in the Landau 
gauge. For a general discussion about lattice Gribov copies see, for example, 1 o and references therein. For 
the SU(2) group, the gauge and ghost propagators versus Gribov copies were studied in The authors 

claim that the gluon propagator is not sensible to Gribov copies in the weak coupling regime^. For the ghost 
propagator, the simulations performed by the first author shows that, close to the continuum, the propagator 
is again not sensible to Gribov copies. In the second study, it is claimed a reduction of 6% for the central 
value of the smaller momenta ghost propagator and a reduction of 4% on the Kugo-Ojima parameter^. The 
SU(2) simulations suggest that the infiuence of Gribov copies is at the level of the simulation statistical 
error. For SU(3) there is no systematic study but it is believed that the Gribov noise is contained within 
the statistical error of the Monte Carlo. 

In this paper we study the pure gauge lattice QCD gluon propagator in the Landau gauge and try to 
understand the role of the Gribov copies. We compute the gluon propagator for 302 configurations, with 
f} = 5.8, for a lattice size of 12**. Our results show that, although being a small effect (less than 6-10%), 
the Gribov copies change the lowest momenta components of the gluon propagator. This effect is illustrated 
fitting the gluon propagator and comparing the results for sets of configurations built from different copies. 

'^Note that, in the strong coupling regime, Cucchieri is able to see differences on the propagator due to Gribov copies. 
^See, also, GHl- 



Gribov copies influence can go from a doubling of the statistical error, to the extreme case of changing the 
functional form of the propagator. 

The paper is organised as follows. Section II sets the field definitions and notation used in this work. 
In section III, the Landau gauge is discussed, both on the continuum and in lattice QCD. Moreover, the 
algorithm used here is sketched. In section IV, the results for the role of Gribov copies in the gluon propagator 
are reported. Finally, in section V our results are discussed. 



2 Field Definitions and Notation 



In the lattice formulation of QCD, the gluon fields are replaced by the links 



(2) 



where are unit vectors along fj, direction. QCD is a gauge theory, therefore the fields related by gauge 
transformations 

Uf^ix) g{x) U^{x) g\x + ae^) , g G 5C/(3) , (3) 

are physically equivalent. The set of links related by gauge transformations to C/^(a;) is the orbit of J7^(a:). 
The gluon field associated to a gauge configuration is given by 



A^{x + ae^/2) = — U^{x) - Ul{x) 
Zigo L 



1 



Qigo 



-Tr 



U,ix) - UUx) 



(4) 



up to corrections of order a^. 

On the lattice, due to the periodic boundary conditions, the discrete momenta available are 

27rn,,, 



9m 



Tin = 0, 1, ... - 1 , 



where is the lattice length over direction /i. Tlic momcntiim space link is 

X 

and the momentum space gluon field 



(5) 



(6) 



U,{q) - UU-q) 



-Tr 



U,{q) - Uli-q) 



(7) 



The gluon propagator is the gluon two point correlation function. The dimensionless lattice two point 
function is 

{Al{q) Al{q')) = D^Uq) V 6{q + q') . (8) 
On the continuum, the momentum space propagator in the Landau gauge is given by 



DfM = S'^" {S,. - ^) D{q') . 



(9) 



Assuming that the deviations from the continuum are negligable, the lattice scalar function D(q^) can be 
computed directly from (O as follows 

Dil') - (jv2„i)(^^,^„i)y E(Tr[^,(g)A,M)]), g^O, (10) 

and 

A'' 

where 

,„^H,n(f), (12, 

Nc = 3 is the dimension of the group, ~ 4 the number of spacetime dimensions and V is the lattice 
volume. 



3 The Landau Gauge 

3.1 The Continuum Landau Gauge 

On the continuum, the Landau gauge is defined by 

d^A^ = 0. (13) 

This condition defines the hyperplane of transverse configurations 

r = {A: d-A = 0} . (14) 

It is well known 0] that T includes more than one configuration from each gauge orbit. In order to try to 
solve the problem of the nonperturbative gauge fixing, Gribov suggested the use of additional conditions, 
namely the restriction of physical configurational space to the region 

n EE {A: d-A^O, M[A] > 0} c T , (15) 

where M[A] = —V • D[A] is the Faddcev-Popov operator. However, is not free of Gribov copies and does 
not provide a proper definition of physical configurations. 

A suitable definition of the physical configurational space is given by the fundamental modular region 
A C ri, the set of the absolute minima of the fmictional 

FA[g] = Id'xY, Tr [Al{x)Al{x)\ . (16) 

The fundamental modular region A is a convex manifold |2L)| and each gauge orbit intersects the interior of 
A only once |21[I22 |. i.e. its interior consists of non-degenerate absolute minima. On the boundary 9A there 
are degenerate absolute minima, i.e. different boundary points are Gribov copies of each other |22[ 1231 
The interior of A, the region of absolute minima of (|16|l . identifies a region free of Gribov copies. 



3.2 The Lattice Landau Gauge 

On the lattice, the situation is similar to the continuum theory |25[ 1261 ITf] , The interior of A consists of 
non-degenerate absolute minima of the lattice version of (|16ll and Gribov copies can occur at the boundary 
9A. However, for a finite lattice, the boundary i9A, where degenerate minima may occur, has zero measure 
and the presence of these minima can be ignored |26j . 

On the lattice, the Landau gauge is defined by maximising the functional 

Fu[g] - CFY,^G{Tr[g{x)U^{x)g\x + fi)]} (17) 

where 

is a normalization constant. Let C/^ be the configuration that maximises -FViff] on a given gauge orbit. For 
configurations near on its gauge orbit, we have 

+ =i^c;[l] + ^^zc^'^(a;)Tr[ {U ^{x) - U ^,{x - (i)) - 

X,fl 

(Ulix) - Ulix ^ f,))] , (19) 

where A° are the Gell-Mann matrices. By definition, C/^ is a stationary point of F, therefore 

dF iCp 



du'^ix) 4 



■^Tr[ X'^iU^ix) - U^ix^fi)) - 

X^{Ulix) - Ulix-fi))] = 0. (20) 
In terms of the gluon field, this condition reads 

^Tr[ A"(A^(a; + a/i/2)- A^(x-a/i/2)) ] + 0{a^) = 0, (21) 

or 

J2d,A;ix) + 0{a) = 0, (22) 

i.e. (|20|l is the lattice equivalent of the continuum Landau gauge condition. The lattice Faddeev-Popov 
operator M{U) is given by the second derivative of (|17|l . 

Similarly to the continuum theory, on the lattice one defines the region of stationary points of (|17|l 

T = {U : d-A{U) ^ 0}, (23) 
the Gribov's region of the maxima of (|17|l . 

n = {U : d- A{U) = and M{U) > 0} (24) 

and the fundamental modular region A defined as the set of the absolute maxima of H17|l . 

A proper definition of the lattice Landau gauge chooses from each gauge orbit, the configuration be- 
longing to the interior of A. 



3.3 Gauge Fixing Algorithm 



On the lattice, gauge fixing is implemented by maximizing Fij[g]. In this work, the gauge fixing algorithm 
used is a Fourier accelerated steepest descent method (SD) as defined in [2H|. In each iteration, the algorithm 
chooses 



g{x) = exp 

where 



(25) 



A_,(C/,,(x)) = U^,{x - ae,) - U^,{x) , (26) 

are the eigenvalues of (—5^), a is the lattice spacing and F represents a fast Fourier transform (FFT). 
For the parameter a we use the value 0.08 |28| . For numerical purposes, it is enough to expand to first order 
the exponential in H25(l . followed by a reunitarization of g{x). 

On the gauge fixing process, the quality of the gauge fixing is measured by 



VN, 

where 



^Tr[A(x)At(x)] (27) 



A(a;) ^ [U^{x-ae^) - Ul{x) - h.c. - trace] (28) 



is the lattice version of d^Afj_ — 0. 



4 The Gluon Propagator 

In this work only pure gauge quenched configurations are considered. The Wilson action configurations were 
generated with version 6 of MILC code ||29j . 

The function Fjj has many maxima - see, for example, j30j . In each gauge orbit, the different maxima 
are different configurations and, therefore, the gluon propagator changes according to the chosen set of 
maxima. In order to study such dependence, 302 gauge configurations were generated for a 12'' lattice and 
for /3 = 5.8, using a combined update of 4 over-relaxed and 5 quasi-heat bath Cabbibo-Mariani updates, 
with a separation between configurations of 3000 combined updates. To each gauge configuration, 500 SD 
gauge fixings, starting from different randomly chosen points, were performed requiring 

^ - ^ETr[A(x)At(a:)] = < 10-^^ (29) 

X X 

From these 500 SD, on each gauge orbit, we keep the gauge configurations associated to the largest maximum 
of Fjj (named MAX in the following), the smallest maximum of Fjj (named MIN) and three random values of 
Fjj (RNDl, RND2, RND3), generated starting the gauge fixing process by choosing always the same random 
g{x) matrices. A further gauge fixing (named ID), starting the gauge fixing process by setting all g{x) = 1, 
was performed to all gauge configurations. Another gauge fixing (named RND), starting the gauge fixing 
process by choosing always the same random g(x) matrices, was performed to all configurations. 
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Table 1: Bare gluon scalar function. The numbers in parentheses are the statistical errors, computed using 
the jackknife procedure, on the last digits of D{q^). 



4.1 Bare Gluon Propagator 

The scalar function D{q^), computed according to equations pO|) and pi|l . after averaging over equivalent 
momenta"^, is shown in figure^as function of q and as function of q. The figures include D{q^) as function of 
momenta of type (g, 0, 0, 0), (q, g, 0, 0), (g, q, g, 0) and (g, q, g, q) for all available g in our lattice. The figures 
for D(q^) are reported in tabled From now on, unless stated clearly, we will consider only the data refering 
to D{q^) as function of q. 

FigureQland table^show that, for the gluon propagator, the effect of Gribov copies is small and visible 
for the smallest momenta. Indeed, comparing the different gluon propagators to the MAX propagator, it 
comes that, within one standard deviation, the ID propagator agrees with the D{q^) MAX for almost all the 



^For example, for each gauge configuration the quoted value for momenta (1, 0, 0, 0) is the average over (1, 0, 0, 0), (0, 1, 0, 0), 
(0, 0, 1, 0) and (0, 0, 0, 1) values. Similarly, for (1, 1, 0, 0) a Z4 average is performed, etc. 



momenta considered. The exception being D{q^) for the momenta associated to = (5, 5, 5, 0), compatible 
with the MAX value only within two standard deviations. Note that only the ID and MAX values agree for 
the infrared regime. The RND propagator agrees, within one standard deviation, with the MAX propagator 
for all momenta but the zero momenta. The zero momenta RND propagator agrees with the MAX D{Q) 
only within two standard deviations. The strongest deviation from the MAX propagator occurs when D{q^) 
is computed using the smallest of the maximum of Fjj. The MIN propagator agrees, within one standard 
deviation, with MAX for momenta |g| > 1.7320 for momenta of type (g, 0,0, 0), \q\ > 2.000 for {q,q, 0,0), 
\q\ > 3.3461 for momenta {q,q,q,0) and \q\ > 3.4641 for (q,q,q,q) momenta. For smaller momenta the 
differences between the D{q^) values can achieve six standard deviations. Indeed, the agreement between 
the MIN and MAX values quoted in the table are: six standard deviations for = (1, 1, 0, 0); four standard 
deviations for = (1,1,1,0) and (2,0,0,0); three standard deviations for — (0,0,0,0), (1,0,0,0) and 
(1, 1, 1, 1); two standard deviations for = (2, 2, 0, 0), (3, 0, 0, 0), (3, 3, 3, 3) and (4, 4, 4, 0). The lattice data 
shows clearly that Gribov copies change the low momenta (q < 1.7320) components of the gluon propagator. 



For zero momentum, the largest propagator occurs when the configurations are gauge fixed to the 
fundamental modular region. The absolute difference between the MIN, RND and ID to the MAX zero 
momenta propagator central values are 6.7%, 2.6% and 1%, respectively. These numbers can be read as an 
order of magnitude of the maximal change on the gluon propagator due to Gribov copies. For the other 
momenta, it is not always true that the largest value of D{q^) is associated to the MAX propagator. This 
can be seen in figure |21 

Figure 121 suggests that the ratio between the propagators to the MAX propagator is a function of q, that 
converges to one for the larger momenta. Moreover, the figure shows clearly that the MIN propagator is 
different from the MAX propagator for momenta smaller than q ^ 1.7. From figure |2] one can quantify again 
the change on the gluon propagator due to Gribov copies. For the MIN propagator, the effect of Gribov 
copies is, at most, a factor of 5-10%. For the RND and ID propagators, the effect of Gribov copies is not so 
dramatic (a factor smaller than 5%). 



Figure |2 could suggest that the ratio between the propagators to the MAX propagator would be a 
constant factor. To test this hypothesis, in figure |3| the propagators are ploted after rescaling the different 
gluon fields to reproduce the central value of the zero momenta MAX scalar function. As seen in figure |31 
the propagators differ by more than one standard deviation for certain momenta. If, instead of rescaling 
the gluon field to reproduce the zero momenta MAX scalar function, the matching is done, for example, for 
= (6, 0, 0, 0), then ^(O) = 36.68 ± 0.46, 36.53 ± 0.50, 35.84 ± 0.47 and 34.70 ± 0.48 for the MAX, ID, 
RND and MIN propagators, respectively. The MAX, ID and RND D{0) are compatible within one standard 
deviation. The MIN D{0) agrees with MAX value only within three standard deviations. In order to further 
test the hypothesis under discussion, D/Dmax was fitted to a constant. No correlations were considered in 
the fits. The /d.o.f. for these fits are reported in figure^J Although, in general, the values of the /d.o.f. 
decrease as one excludes more lower momenta, they are always too high to conclude that the ID, RND and 
MIN propagators differ, from the MAX propagator, by a multiplicative factor. In particular, the difference 



between MIN, ID and MAX propagators is clearly not a constant. The RND/MAX ratio is compatible with 
a constant for the largest momenta considered. 

In conclusion, the analysis of the raw data for the bare gluon propagator suggests that the effect of 
Gribov copies is small, but observable (clearly, less than a 10% factor) and is stronger for smaller momenta. 
Moreover, Gribov copies have almost no effect on the high momentum components of the gluon propagator. 
The data reported in tabled shows that the effect of Gribov copies can be overcame if one multiplies the 
statistical errors by a factor of 2 to 3 for the smaller momenta {aq < 1.73). This doubling of the statistical 
error can be either, a general property associated to the effect of Gribov copies, or a result due to the limited 
statistics used here. Note that in the SU(2) study of 16' , the number of configurations used for the larger 
lattices (12**, 16^) and for the larger f3 {— 2.7) was about half or less than half of the configurations used in 
our simulation. The investigation of D{q^) / Dy^jj^y^{q^) shows that the propagators associated to the Gribov 
copies named as ID, RND and MIN differ from the MAX propagator by more than a constant factor. 



4.2 Gribov copies and gluon propagator models 



In the previous section, it was argued that the ID, RND and MIN propagators do not differ from the MAX 
propagator by a constant factor. The question we would like to investigate now being: is it possible to 
quantify the differences, due to Gribov copies, when modelling the gluon propagator? To try to answer this 
question, we will study the fit of D(q^) to a functional form. 

In |31| a number of gluon propagator models were studied. Our simulation access a limited range of 
momenta and, certainly, finite space and/or finite volume effects are no negligeable. Instead of performing a 
detailed study of several functional forms, we chose to investigate the model which, according to Lcinweber 
et al describes better the lattice data. 

Let us assume that the scalar function is given by 



where 



D{q^) 



L{q\M^) 



[q^ + AP) 



l+a 



L{q\M^) 
q^ + AP 



ln[(g2 + Af2)(g-2 + M-2)] 



(30) 



(31) 



is an infrared-regulated version of the one-loop logarithm correction to the gluon propagator and, for pure 
gauge theories, = 13/22. 

According to the results of the previous section, Gribov copies seem to change the gluon propagator 
for the low energy momenta. Therefore, to measure such an effect we will consider three different types of 
uncorrelated fits. A fit to the highest momenta (UV-fit) using the following functional form 



Diq') 




(32) 



a one-loop corrected perturbative gluon propagator. A fit to the lowest momenta (IR-fit), assuming that 



D{q' 



(g2 + M2) 



l+a 



(33) 





Z A x^ld.o.j. 


MAX 


1.4731}} 0.80761^^ 0.08 


ID 


1.4578j:^^ 0.818lt^^ 0.48 


RND 


1.4620t;^^ 0.8167t^? 0.44 


MIN 


1.4243+^^ 0.8465l[!^ 2.78 



Table 2: Fits of the higher momenta of type = (n, 0,0, 0) to the one-loop corrected perturbative gluon 
propagator (|32|l . The fitting range goes from n = 3 to n = 6. For larger fitting ranges, the /d.o.f. becomes 
too large (> 18). It is possible to fit the data using a smaller fitting range (n = 4 to n = 6). However, 
we do not report the figures because such a fit would have only one degree of freedom. Statistical errors 
were computed using the bootstrap method. The quoted errors correspond to a 68% confidence limit of the 
distributions obtained from 5000 bootstrap samples. 





= (2,2,0,0) 
Z A xVd.o.f. 


= (3,3,0,0) 
Z A xVd.o.f. 


MAX 
ID 
RND 
MIN 


1. 78461^^ 0.7208t^'j 1.44 
1.7867t^4 0.7184t^^ 0.44 
1.7776t^g 0.7285^33 0.82 
1.7617t^^ 0.7413^^? 0.52 


1.86lt}^ 0.659t}^ 0.14 
1.823t}^ 0.688t}^ 0.15 
1.829t}? 0.686j:}i} 0.27 
1.804t}^ 0.706t}^ 0.72 



Table 3: Fits of the higher momenta of type = (n,n, 0,0) to the one-loop corrected perturbative gluon 
propagator H32(l . The fitting range goes from n = 2 or 3 up to n = 6. Statistical errors were computed using 
the bootstrap method. The quoted errors correspond to a 68% confidence limit of the distributions obtained 
from 5000 bootstrap samples. 

and a fit of H30|l to all lattice data. 

In order to compare our results with we take their central values for at /3 = 6.0 and /3 — 6.2 and 
scale a to /? = 5.8 using the results of two-loop calculations. This procedure gives, respectively, — 1.463 
GeV and a^^ = 1.590 GeV. The average of the two values being a^^ = 1.53 ± 0.06 GeV (a = 0.13 fm). 



Figure |31 shows q^D{q^) as function of q for all sets of gauge fixed configurations. The results for the 
different momenta shows that, in our simulation, the finite space/volume effects are not negligable - an effect 
of the order of 10% from ((7,0,0,0) to the other types of momenta. Since the different types of momenta 
have different finite space/volume effects, we will not include different types of momenta in the fits. The 
exception being the IR fits. 

The fits of the highest momenta to the asymptotic form 1)3211 are reported in tables|21to[Slfor all types 
of momenta. The first point to remark is that the gluon propagator scales perturbatively for aq > V2 for 
momenta associated to = (n, 0,0,0) and ~ (n, n,0, 0), for aq > 2.450 for — {n,n,n,0) momenta 
and for aq > 1.035 for = {n,n,n,n) momenta; i.e. the asymptotic form describes quite well the lattice 
data for sufficiently large momenta. Perturbative scaling starts at momenta q ^ 1.6 — 3.7 GeV, a value 
compatible with the figure quoted in 31 , 2.7 GeV. 

In what concerns the effect of Gribov copies at high momenta, the results given in tables 12 to ^ show 





Z A x^ld.o.j. 


MAX 
ID 

RND 
MIN 


2.100t^^ 0.534t^^ 0.78 
2.165i^5 0.498t^^ 1.29 
2.1501:3^ 0.499^27 0.09 
2. 09211^^ 0.555^27 0.80 



Table 4: Fits of the higher momenta of type = (n, ri, n, 0) to the one-loop corrected perturbative gluon 
propagator The fitting range goes from 3 up to n = 6. For larger fitting ranges, the /d.o.f. becomes 
too large (> 2). It is possible to fit the data using a smaller fitting range (n = 4 to n = 6). However, 
we do not report the figures because such a fit would have only one degree of freedom. Statistical errors 
were computed using the bootstrap method. The quoted errors correspond to a 68% confidence limit of the 
distributions obtained from 5000 bootstrap samples. 





n,, = (1,1,1,1) 
Z A .'^^ 

d.O.J. 


n^-(2,2,2,2) 
Z A 

d.O.J. 


MAX 
ID 

RND 
MIN 


2.099+^^1 0.5944+*^ 0.39 
2mQt\i 0.60431^? 0.52 
2.0756^^? 0.6122l4^ 1.12 
2.0146+^2 0.65161:^'^ 1.20 


2.112l^( 0.584+^^ 0.51 
2.102^27 0.584t^^ 0.64 
2.064t^^ 0.621^20 1-48 
2.09lti^j^ 0.590tg 1.05 



Table 5: Fits of the higher momenta of type — {n,n,n,n) to the one-loop corrected perturbative gluon 
propagator H32I) . The fitting range goes from n = 2 or 3 up to n = 6. Statistical errors were computed using 
the bootstrap method. The quoted errors correspond to a 68% confidence limit of the distributions obtained 
from 5000 bootstrap samples. 



that, for the same data and fitting range, the MAX, ID and RND values are compatible within one standard 
deviation. For momenta associated to = (n,n,n,7i) and for the largest fitting range**, the Z and A 
values are compatible within two standard deviations. On the other side, the MIN fitted parameters are 
not compatible with the MAX figures; the exception being the fit to (g, q, q, 0) momenta and the fit to the 
smallest fitting range reported in tabled 

In what concerns the stability of results, in general, the fitted parameters are not stable against a 
change in the fitting range. Probably, this is due to the limited number of different momenta available in 
the simulation. If one compares the results of the larger fitting ranges where Z and A are compatible within 
one standard deviation for the different types of momenta, it comes that Z increases and A decreases as one 
goes from rt^ = (rt, 0, 0, 0) to — (ti, n, n, n) by a factor of ~ 1.4. Such a large correction is an indication 
of important finite space effects - remember that the lattice spacing is 0.13 fm. If one compares our values 
for A with those reported in [2], the numbers given in tablesEltoElare, typically, larger than those reported 
by Leinweber et al. 

The discussion of the IR properties of the gluon propagator requires data for small momenta. In our 
simulation one has only a limited access to the infrared regime of QCD. This is a serious limitation to a 
proper investigation of the low energy gluon propagator. Nevertheless, we have tried to find the combination 
of the smaller momenta which is well reproduced by (|^ . Unfortunately, to achieve such a goal, we had 



*In physical units, the fitting range includes momenta from 1.6 GeV up to 6.1 GeV. 





A 


M 


a 


xVd.o.f. 


MAX 
ID 

RND 
MIN 


17.72+26 
18.19^29 

18.78^2? 
22.8lt^3 


0.6947lg 
0.7076t^^ 
0.7237t^^ 
0.8189t^g 


1.278i^3 
1.312+2^ 

1.363t24 
1.6751^7 


0.048 
0.031 
1.032 
1.561 



Table 6: The infrared propagator. Statistical errors were computed using the bootstrap method. The quoted 
errors correspond to a 68% confidence limit of the distributions obtained from 5000 bootstrap samples. 



to combine different types of momenta. Below, we will show the results of such investigation. The reader 
should be aware of the physical meaning, or lack of meaning, of the numbers reported here. We would like 
to remember that our main goal is to see if there are differences, on the propagators, due to the choice of 
Gribov copies. 

The set of momenta associated^ to = (0, 0, 0, 0), (1, 0, 0, 0), (1, 1, 0, 0) and (2, 0, 0, 0) is well described 
by the model function H33|l . The fitted parameters are reported in table El for the different propagators. 
Although, the lattice data is well described by (|33|l . not all fitted parameters are compatible within one 
standard deviation. Indeed, the MIN propagator values are not compatible with any of the other propagators. 
The MAX and ID propagators all have the same A parameter. The A from the RND fit is, within two 
standard deviations, compatible with the MAX figures. In what concerns the gluon mass M, the MAX and 
ID values are compatible within one standard deviation but MAX and RND are compatible within three 
standard deviations. For the parameter a, the MAX and ID values are compatible within one standard 
deviation but MAX and RND are compatible within two standard deviations. Note that the gluon mass M 
computed from the IR regime of QCD is not compatible, within one standard deviation, with the values of 
A from the UV regime - see tables |21 to [S] The values of M and a for MAX are the smallest figures in table 
El Prom these fittings, one can quantify the effect due to Gribov copies as a two to three sigma effect on the 
parameters. 

Finally, let us discuss the fittings of (|30|) to all lattice data. The results of the fittings are reported 
in table for momenta — (n, 0,0, 0), in table [S] for — 0,0), in table El for = (n, n,n, 0) 

and in table 11(11 for = {n,n,n,n) momenta. The x^/d.o.f. shows that, in general, the lattice data is 
well described by 1)30(1 . The exceptions are the fits to the MIN data, momenta = (n, 0,0,0), and ID 
propagator, momenta — n,0). For these two cases the x^/d-o.f. is quite large, meaning that the 

lattice data is not described by 130|) . 

To identify the effect of Gribov copies the different fittings are compared for the same type of momenta. 
The data on tables[3toEIlshows that, for all types of momenta, the fitted parameters for the MIN propagator 
are not compatible with the corresponding parameters for the MAX propagator. For momenta associated 
to = (n, 0, 0, 0), the ID and MAX propagators parameters are compatible within one standard deviation. 
The RND and MAX Z values are compatible, within the same level of precision, the a and M values are 
compatible within 2cr and A is compatible within three standard deviations. For — (n, n, 0,0) momenta, 
ID and RND parameters are compatible with the MAX values only within two standard deviations. The 
exception being the a from RND propagator, which agrees with the MAX figures within la. For = 



= 0, 0.52, 0.73 and 1 or, in physical units, q = 0, 0.80, 1.12 and 1.53 GeV, respectively. Note that the number of degrees 
of freedom for this fit is one. 



Z A M a x^/d.o.f\ 

MAX] 1.58lt};^ 12.631^9 0.737^^^ 1.982+56 

ID 1.564tJi 13.04t^? 0.748i" 2.004lg 0.39 

RND 1.580tio 13.82t^? 0.780tii 2.1341^^ o.74 

MIN 1.559^*;^ 15.691^2 0.84lt^^ 2.320±^g 3.45 



Table 7: Fits to all lattice data for momenta associated to rt^ = (n, 0,0,0) to the fmictional form H30|l. 
Statistical errors were computed using the bootstrap method. The quoted errors correspond to a 68% 
confidence limit of the distributions obtained from 5000 bootstrap samples. 





Z 


A 


M 


a 


xVd.o.f. 


MAX 
ID 

RND 

MIN 


1.85651^1 
1.8478t^^ 
1.8430t|8 
1.8055li^^ 


10.46lf^ 
11.021}^ 
10.96li^ 
13.94j:^^ 


0.7283^'^^ 
0.7493l:^g 
0. 75241;^^ 
0.8569t^8 


1.990+^1 

2.061^29 
2.046+^° 

2.389^31^ 


1.15 
0.12 
0.45 
0.33 



Table 8: Fits to all lattice data for momenta associated to — (n, n,0, 0), with n from to 6, to the 
functional form Ip^Hjl . Statistical errors were computed using the bootstrap method. The quoted errors 
correspond to a 68% confidence limit of the distributions obtained from 5000 bootstrap samples. 



(n, n,ri, 0), RND and MAX values are compatible within two standard deviations. For the ID parameters, 
the Z value is, within two standard deviation, compatible with the MAX value and all remaining parameters 
are compatible within la. For = (n,n,n,n), the MAX, RND and ID fitted parameters are compatible 
within one standard deviation; the Z for the ID and MAX are compatible within 2a. Note that, in general, 
the MAX propagator has the larger Z value and the smallest A, M and a. Again, like in the IR fits one can 
quantify the effect due to Gribov copies as a two to three sigma effect. From the fittings it is not possible to 
establish, clearly, which parameters are less sensible to Gribov copies. Note that the fits to n^^ = (n, n, n, n) 
momenta, although having large statistical errors and with the exception of the MIN propagator, they do 
not distinguish the Gribov copies. 

From the analysis of tables 171 to 1101 one can check which parameters are robust against change of fitting 
momenta. Indeed, looking only to the fundamental modular region propagators, it comes that the overall 
normalization parameter Z is a function of the type of momenta considered. At lu, the different Z values are 
not compatible with each other. For the same level of precision, the — (rt, n, 0,0), ti^ — (n, 7i,rt, 0) and 
= (n, n, n, n) fitted parameters which measures the relative normalization of the infrared to ultraviolet 
propagator components. A, are compatible with each other. On the other hand, M and a parameters are 
robust against change of momenta. All four values reported in the tables are compatible within one standard 
deviation. Our results for a and M are 





(g, 0,0,0) 


(9,9,0,0) 




{q,q,q,q) 


a 

M (MeV) 


1.982l^« 
1128li^±44 


1.990+i^i 
11141:^^ _^ 44 


2.018+:!^ 
1128ti4±44 


2.1511^ 
11571^^ ±45 



where the second error in M represents the error in the lattice spacing. Curiously, these values are compatible, 
within one standard deviation with the values quoted in Jl,, namely a = 2.21q 21o'.3 ^ = (1020 ± 100 ± 



MAXI 
ID 

RND 

MIN 



Z 



A 



M 



a x^/d.o.f. 



1.94101!^^ 10.26 
1.95351!^? 
1.9289^4^ 
1.90901-^^ 



IT 

-27 

10.40t25 

12.83i26 



0.7371 
0.7484 
0.773 



0.8460i;|5 



-91 
+97 

-89 
+ 10 



2.018^39 

2.07lt^i 
2.144^41 

2 .355 



1.71 
2.46 
1.60 
0.68 



Table 9: Fits to all lattice data for momenta associated to — {n,n,n,0), with n from to 6, to the 
functional form H30() . Statistical errors were computed using the bootstrap method. The quoted errors 
correspond to a 68% confidence limit of the distributions obtained from 5000 bootstrap samples. 



MAXi 
ID 

RND 
MIN 



A 



M 



a X /d.o.f 



-13 
-11 
^-13 

2.001^}^ 

-70 
-87 



2.018: 
1.993^ 



1.9239: 



10.38+8^ 

11.07j:^^ 
10.42t^2 

+74 
-69 



0.756: 



-28 

0.780124 
0.764124 



14.30 



0.896: 



-21 
-20 



2.15 
2.22 
2.117 
2.535 



+TD- 

-12 
+ 10 
-11 
+98 
-95 
+89 
-86 



0.48 
0.44 
1.12 
0.57 



Table 10: Fits to all lattice data for momenta associated to — {n,n,n,n), with n from to 6, to the 
functional form H30(l . Statistical errors were computed using the bootstrap method. The quoted errors 
correspond to a 68% confidence limit of the distributions obtained from 5000 bootstrap samples. 



25) MeV. Note that the values for M and a quoted for the fittings to all lattice data are not, in general, 
compatible with the same parameters computed from the IR and UV fits. Probably, this is due to using 
a relatively small lattice that does not enable a clear separation between the low energy and high energy 
regimes of QCD. 



5 Discussion and Conclusions 

In this paper the problem of the Gribov copies in lattice QCD is addressed. To try to understand the role of 
Gribov copies in lattice QCD, the gluon propagator was computed with 302 configurations for a 12^ lattice 
and for (3 = 5.8 using the overrelaxed quasi-heat bath. 

The analysis of the raw data shows that Gribov copies change only the low momenta components of 
the gluon propagator. In our simulation, only for momenta aq < 1.7320 (g < 2.6 GeV) there are significant 
differences between the MIN and MAX propagators. The RND data is not compatible, within one standard 
deviation, with the MAX propagator only for zero momentum. The study performed here shows that, 
typically, the choice of different Gribov copies changes the propagator in such a way that the figures become 
compatible within two-to-three standard deviations. Note, however, that for the patological case of the 
MIN propagator the deviation relative to the MAX propagator, can be as large as six standard deviations 
effects. This result seems to suggest that in the study of the IR limit of the gluon propagator, the statistical 
errors should be multiplied by a factor of two or three in order to take into account possible deviations due to 
Gribov copies. If this is true for the statistical accuracy of our study, this may not hold when larger statistics, 
bigger lattices are considered. That the Gribov copies change essentially the IR limit of the propagator can 
be seen in figure El where D{q^) is plotted against {Fu). 



The properties observed for the raw scalar data are observed when we model the lattice data. A 
difference, due to the choice of Gribov copies, of up to two-to-three standard deviations is seen on the fitted 
parameters. Note that this is observed even if the lowest energy momenta have the largest absolute errors, 
i.e. there contribution to the is not so relevant. To our mind, a deviation of this order of magnitude is, 
probably, a good measure of the influence of the Gribov copies on the gluon propagator. 

In what concerns the gluon propagator itself, the results of our simulation for M and a support the 
results quoted in a previous investigation using much larger lattices Remember that, for the MAX 

propagators, these parameters are robust against a change on the type of fitted momenta. It is interesting, 
that the lattice data supports quite well a gluon propagator which, for large momenta, behaves like a massive 
vector with a mass of the order of the hadronic scale. Remember that a massive gluon propagator, with a 
gluonic mass of the order of 1 GeV, has some phenomenological support |32|. 

P. J. S. acknowledges financial support from the Portuguese FCT. This work was supported by grant 
SFRH/BD/10740/2002. 
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(a) Scalar Function as function of q. 
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(b) Scalar Function as function of q. 



Figure 1: Bare gluon propagator. Statistical errors were computed using the jackknife procedure. 
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Figure 2: D{q^) / DMAxiff') as function of (7 for ID, RND and MIN propagators. Statistical errors were 
computed using the bootstrap method. The quoted errors correspond to a 68% confidence Umit of the 
distributions obtained from 5000 bootstrap samples. 



MAX 
ID 

RND 
MIN 



1^ -L 



(a) D{q^)D-^ji^y^{0) / D{0) as function of q for all momenta. 




(b) D{q^)D^j:^y^{0) / D{0) as function of q for the larger momenta. 



Figure 3: Scaled gluon propagator. Statistical errors were coniputed using the bootstrap method. Statistical 
errors were computed using the jackknife procedure. 
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Figure 4: /d.o.f. for the uncorrelated fits of D{q^) / D MAxio^) to a constant. In the figure, i is the number 
of lower momenta, not considered in the fit. 




Figure 5: q^D{q^) as function of q for MAX, ID, RND and MIN propagators. The points not connected 
by lines refer to (g, 0,0,0) momenta. The points connect by dotted hues refer to (g, (7,0,0) momenta, the 
points connected by dashed Unes to (g, q, q, 0) momenta and the points connected by dashed-dotted hnes to 
(g, q, q, q) momenta. Statistical errors were computed using the jackknife procedure. 
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(a) Zero momentum scalar function. 



(b) Scalar function associated to Up, = (1, 0, 0, 0). 
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(c) Scalar function associated to = (1, 1,0,0). 



(d) Scalar function associated to = (1, 1, 1,0). 



Figure 6: Evolution of D{q^) with F. After ordering the gauge fixed configurations associated with the 
sets MAX, ID, RND, MIN, RNDl, RND2 and RND3 according to the Fu value, the gluon scalar function 
is computed picking always configurations within the same class of values of Fu. Statistical errors were 
computed using the jackknife procedure. 



